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ABSTRACT 

In the unification scheme, narrow lined (type 2) active galactic nuclei (AGN) are intrinsically 
similar to broad lined (type 1) AGN with the exception that the line of sight to the broad 
emission line region and accretion disk is blocked by a dusty torus. The fraction of type 1 AGN 
measures the average covering factor of the torus. In this paper, we explore the dependence 
of this fraction on nuclear properties for a sample of low redshift (z < 0.35) radio strong 
CPi.4GHz 5 s 10 23 W Hz -1 ) AGN selected by matching the spectroscopic catalog of Sloan Digital 
Sky Survey and the radio source catalog of Faint Image of Radio Sky at Twenty cm. After 
correcting for several selection effects, we find that : (1) type 1 fraction f\ keeps at a constant 
of ~ 20% in the [Om] luminosity range of 40.7 < log(L[om]/erg s -1 ) < 43.5. This result is 
significantly different from previous studies, and the difference can be explained by extinction 
correction and different treatment of selection effects. (2) f\ rises with black hole mass from 
~ 20% below 1O 8 M0 to 30% above that. This coincides with the decrease of the fraction of 
highly-inclined disk galaxies with black hole mass, implying a population of Seyfert galaxies 
seen as Type-2 due to galaxy-scale obscuration in disk when the host galaxy type transfer 
from bulge-dominant to disk-dominant. (3) f\ is independent of the Eddington ratio for its 
value between 0.01 and 1; (4) f\ ascends from 15% to 30% in the radio power range of 
23 < log(Pi.4GHz/WHz _1 ) < 24, then remain a constant at ~ 30% up to 10 26 W Hz" 1 . 

Key words: quasars: general — galaxies 



1 INTRODUCTION 

Active Galactic Nuclei (AGN) are traditionally divided into 2 sub- 
classes: type 2 and type 1, according to the absence or presence 
of the broad emission lines. Their different observed properties can 
be explained largely via anisotropic obscuration, likely by a dusty 
torus on scales of several to tens parsecs, of otherwise the same 
type of objects (Antonucci 1993; Tran 1995; Nenkova et al. 2002). 
Because of its large physical scale, Narrow emission Line Region 
(NLR) cannot be (completely) obscured by such torus. Therefore, 
an AGN will appear as a type 2 source, when our line of sight to 
the Broad emission Line Region (BLR) and the accretion disk is 
blocked by the dusty torus. Scattered broad lines have been detected 
in the polarized light in almost half of type-2 AGN, which strongly 
supports the unification scheme (Antonucci & Miller 1985; Miller 
& Goodrich 1990; Moran et al. 2000). 

An important parameter in this model is the opening angle of 
the torus. The fraction of type 1 AGN (hereafter f\) is a measure 
of the average opening angle of torus in such a unification scheme. 
Previous studies have shown a type 2 to type 1 ratio of 3: 1 for local 
Seyfert galaxies (e.g., Maia et al. 2003; c.f, Ho et al. 1997). How- 
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ever, there is no reason that the torus opening angle should be the 
same for every AGN, rather it may depend on the black hole mass, 
accretion rate, luminosity or other intrinsic parameters. Exploring 
these parameters dependence will be an important extension to the 
simple unification scheme, and yield the insight into the origin of 
the dust torus as well. 

The luminosities of narrow emission lines, such as [O m], mid- 
infrared light and hard X-rays have been considered as isotropic 
properties (e.g. Meisenheimer et al.2001; Kauffmann et al. 2003; 
Wang et al. 2006). Thus, we can examine the luminosity depen- 
dence of /i using those luminosities. It was found that f\ increases 
with [O m] luminosity in a large sample of Seyfert galaxies from 
Sloan Digital Sky Survey (SDSS, York et al 2000) Data Release 2 
(DR2) (Simpson 2005; Hao et al. 2005). These results were inter- 
preted in the context of "receding torus" (Lawrence 1 99 1 ; Hill et al. 
1996). Similar results have been reported from statistical studies of 
hard X-ray selected AGN (Ueda et al. 2003; Hasinger et al. 2004; 
Gilli et al. 2007; Hasinger et al. 2008; c.f, Wang et al. 2006). 

Radio loud AGN, including radio quasars and radio galaxies, 
can also be unified in such scheme (see Urry & Padovani 1995 for 
a review). As in the case for Seyfert galaxies, broad permitted lines 
were detected in the polarized light of narrow line radio galaxies 
(e.g., Antonuuci, Hurt & Kinney 1994; Tran, Cohen & Goodrich 
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1995; Young et al.1996; Tran et al. 1998). Barthel (1989) found 
that radio galaxies have larger linear dimensions than that of radio 
quasars in their 3CRR sample, and hypothesized that powerful ra- 
dio galaxies and radio quasars belong to the same population, with 
a cone angle of 40° - 50° to mark the division between the radio 
quasar and galaxy. Within this scheme, by using 172 3CR radio 
sources, Lawrence (1991) found that the type 2 fractions are anti- 
correlated with radio luminosity in the range of L llmll2 > 10 25 
W Hz -1 . Falcke et al. (1995) proposed that for radio-loud sources, 
jets may clear a path through the dust, and cause the obscuration 
along the jet's periphery, and give rise to the anti-correlation be- 
tween obscuration and radio luminosity. These mentioned above 
only focused on the most powerful radio AGN. 

In this paper, we will extend those analysis to radio 
moderately-strong AGN, and examine in more details f t as a func- 
tion of radio power, [O m] luminosity, black hole mass and accre- 
tion rate using the large sample of AGN in the SDSS spectroscopic 
catalogues of galaxies and quasars that have been detected in Faint 
Image of Radio Sources at Twenty cm ( FIRST, White et al.1997). 
A combination of large sky coverage, high completeness to a rela- 
tively deep magnitude for galaxies and quasars, and moderate spec- 
tral resolution makes the sample ideal for such a study. 

Based on the SDSS spectroscopic data set, we culled 711 
type 2 objects and 286 type 1 objects having FIRST luminosity 
^i.4GHz > 10 23 W Hz -1 and z < 0.35. After correcting several se- 
lection effects, we obtained f\ as a function of nuclear parameters. 
The paper is arranged as follows. The sample is described in the 
next section. The selection function is estimated in §3 and the re- 
sults are presented in §4. Finally, we will discuss our result in §5. 
Throughout this paper, we will adopt a concordant cosmology with 
H = 71 km s _1 Mpc- 1 , Cl m = 0.27, and Q,, = 0.73 (Spergel et al. 
2003). 



2 THE SAMPLE OF RADIO STRONG AGN 
2.1 Parent Sample 

Starting with the spectroscopic samples of quasars and galaxies 
in the SDSS data release four (DR4), we construct the low red- 
shift sample of radio detected galaxies and AGN. A redshift cut 
z < 0.35 is applied so that Ha falls in the SDSS spectral cover- 
age. To simplify the estimation of selection effects, we consider 
only the objects targeted as main galaxies (Petrosian magnitude 
r < 17.77 and the target mask as "TARGET.GALAXY") or low 
redshift quasars (psf magnitude i < 19.1 and the target mask 
as "TARGET_QSOs_CAP/SKTRT"), or FIRST counterparts (unre- 
solved objects with FIRST counterparts within 2" and psf magni- 
tude ; < 19.0, masked as "TARGET_QSOs_FIRST_CAP/SKIRT"). 
This low-z galaxy and quasar sample (405,904 SDSS spectra in 
total) are then cross-correlated with the FIRST source catalog 
(Becker et al. 2003) to form the radio detected galaxy and AGN 
sample with a procedure described in Lu et al. (2007). 

We use positional coincidence to select radio point sources, 
and visually inspect the FIRST images for all the candidates with 
extended radio morphology (refer to Lu et al. 2007 for details). 
Briefly, we take 2'.'0 as the cutoff of position offset for SDSS-FIRST 
positional coincidence, which is a trade-off between the complete- 
ness and random contamination. Then, we visually inspected the 
cutouts of FIRST image to select the extended lobe(s) apart from 
SDSS nucleus, using a degraded image of 3CR radio sources as the 
reference for physical association. For z > 0.073, we checked the 



FIRST cutouts in 6' x 6', which corresponds to a physical size of 
500 kpc at the redshift z=0.073, because only a small number of 
radio sources have linear sizes larger than this. For z < 0.073, we 
check the cutout with an angular size corresponding to their linear 
sizes of 500 kpc at the sources' redshift. As the outcome of the sec- 
ond stage, we obtain 20,334 low redshift (z < 0.35) objects from 
SDSS-FIRST matching process. 

Next, we apply a radio luminosity cutoff to the sample. Ac- 
cording to Yun et al. (2001) and Hopkins et al. (2003), star forma- 
tion (SF) galaxies rarely have radio powers larger than Pi aghz = 
10 23 W Hz" 1 . Best et al. (2005) found that the radio luminosity 
function at 1 .4 GHz of SF galaxies intersect with radio loud AGN 
at ~ 10 23 W Hz -1 . Beyond this value, the luminosity function of 
SF galaxies drops dramatically, while that of radio AGN decreases 
mildly. Although at ~ 10 24 W Hz -1 , SF galaxies account for nearly 
10% of the population, they will be rejected on the emission line- 
ratio diagrams (hereafter, BPT diagram; Baldwin, Phillips & Ter- 
levich 1981). So we will limit our analysis to the sources with radio 
luminosities above 10 23 W Hz -1 . 

A ^-correction to the radio luminosity is applied by assuming 
a radio spectral index a = 0.5 (/,, oc v~"). After retaining only the 
highest S/N spectrum for the duplicated observations, we obtain 
7,810 SDSS radio loud sources, 276 from SDSS quasar sample and 
7,534 from SDSS galaxy sample, at redshift z < 0.35 and with 
/'i.4GHz > 10 23 WHz-'. 

Our criteria are similar to Best et al. (2005), who cross- 
correlated SDSS-FIRST "compact" sources with 3" and selected 
the "extended" radio counterparts within 30" by assuming that ra- 
dio sources are "double lobe" or "core-lobe", but extend to a fainter 
limit and also include larger radio sources. These authors extracted 
2,215 radio AGN and 497 SF galaxies brighter than 5 mJy (corre- 
sponding to 10 23 W Hz" 1 at z -0.1) from 212,000 z < 0.3 SDSS 
DR2 "main galaxies" spectroscopic targets, and their AGN classi- 
fication was based on the £)„(4000) versus Li.4ghz/AZ» diagnostic 
plane. 

2.2 Continuum Subtraction and Emission-Line 
Measurements 

In order to classify and obtain the intrinsic properties of the radio 
galaxies and AGN, precise measurements of emission line fluxes 
are necessary. We use the measured parameters from the Value- 
added Extra-GAlactic Catalog developed and maintained by Center 
for Astrophysics, University of Science and Technology of China 
(USTC-VEGAC; X.-B. Dong et al. in preparation). We will de- 
scribe briefly the steps relevant here and leave details to the re- 
ferred paper. We correct the SDSS spectra for the Galactic extinc- 
tion (Schlegel, Finkbeiner & Davis 1998) using the extinction curve 
of Fitzpatrick (1999). The spectra are then brought into their rest 
frame using the redshift provided by SDSS pipeline. The contin- 
uum subtraction and emission line measurements are done sepa- 
rately for the different type of objects as follows: 

The continuum subtraction is done according to their relative 
contribution of star-light and AGN continuum. For the narrow line 
objects, the continuum is dominated by star-light, thus, can be mod- 
eled with Independent Component templates (IC templates) fol- 
lowing the procedures described in Lu et al. (2006). Q In brief, we 

1 IC templates were derived from Independent Component Analysis 
(ICA), developed by Lu et al. (2006). Using this technique, Lu et al. (2006) 
compressed the synthetic galaxy spectral library to six nonnegative inde- 
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fit galaxy spectra with the templates derived by applying Essem- 
bling Learning for Independent Component Analysis (EL-ICA) to 
the simple stellar population library (Bruzual & Chariot 2003). The 
templates were then broadened and shifted to match the stellar ve- 
locity dispersion of the galaxy. In this way, stellar absorption lines 
are reasonably modeled to ensure the reliable measurement of weak 
emission lines. At the same time, stellar velocity dispersion, a cor- 
rection to the redshift as well as an average internal extinction to 
the stellar light is obtained. 

For nucleus dominated type 1 AGN where Fe n multiplets and 
other broad emission lines are highly blended, we fit simultane- 
ously the nuclear continuum, the Fen multiplets and emission lines 
(see Dong et al. 2008 for details). The nuclear continuum is approx- 
imated by a broken power-law. Fe n emission, both broad and nar- 
row, is modeled using the Fe n templates provided by Veron-Cetty 
et al. (2004). For those Seyfert Is with significant contribution of 
starlight as measured by the equivalent widths (EWs) of the Ca n k 
/13934 or high order Balmer absorption lines or Nai ,U5890, 5896, 
a starlight model is also included using the 6 IC templates as de- 
scribed above. The decomposition of host-galaxy starlight, nuclear 
continuum and Fen emission were carried out following the proce- 
dure as described in detail in Zhou et al. (2006). 

After subtracting the continuum, we fit emission lines with 
multi-gaussian model using the code described in detail in Dong 
et al. (2005, 2008). Briefly, each line is fitted with one or more 
Gaussians as statistically justified (mostly with 1-2 Gaussians); 
the line parameters are determined by minimizing^ 2 . The [Om] 
/1/14959, 5007 doublet are assumed to have the same profiles and 
redshifts; likewise, [Nn] /U6548, 6583 and [Sn] /U6716,6731 
doublet are constrained in the same way. Furthermore, the flux ra- 
tios of [O m] doublet and [N n] doublet are fixed to the theoretical 
values. Usually, Ha and [Nn] doublets are highly blended and thus 
hard to be isolated; in such cases, we fit them assuming they have 
the same profile as [S n] doublets, which is empirically justified 
(e.g., Filippenko & Sargent 1988; Ho et al. 1997; Zhou et al. 2006). 
For the possible broad Ha and H/? lines, we use multiple Gaussians 
to fit them, as many as they could be statistically justified. If a broad 
emission line is detected with S/N > 5, we regard it as genuine. If 
the broad H/3 line is too weak to achieve a reliable fit, we then re-fit 
it assuming that it has the same profile and redshift as the broad Ha 
line. 

Examples of emission line modeling in Hff+[Nn] and 
HyS+[Om] regions are illustrated in Figure [Tj The two upper pan- 
els shown the objects which was classified by SDSS as galaxies, 
and the two lower ones shown the SDSS classified quasars. The to- 
tal emission line flux is estimated by adding different components, 
while the line width is measured in the combined model. 



2.3 Optical Spectral Classification 

Based on the emission line parameters measured in the last section, 
we classify optical spectra into broad lined and narrow lined AGN, 
and SF galaxies, respectively. Detail criteria for each class are as 
follows: 

Broad line AGN are culled according to the following crite- 
ria (cf. Zhang et al. 2008, X.-B. Dong et al. in preparation): (1) 
The broad Ho- component has been detected with S/N > 5; (2) The 
height of broad Ha is more than twice of root mean square (RMS) 

pendent components (ICs), which were proved to be good templates for 
modeling most normal galaxy spectra. 




spSpcc — 51941— 0290 — 021. fit. gz 



Figure 1. Emission line fitting in the Ha and HyS regions. The black lines 
represent the observed "continuum free" emission lines and noise. The 
blue lines represent for the fitted broad components and the green lines for 
broad+narrow fitting. 
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Figure 2. The BPT diagram for narrow emission line galaxies with radio 
luminosity Pi.4ghz > 10 23 W Hz~' and z ^ 0.35. The solid lines sepa- 
rate Seyfert from SF galaxies, and the dashed lines divide the "composite" 
galaxies from SF galaxies. The blue line represents for Seyfert/LINER sep- 
aration. 

of the continuum-subtracted spectrum in the neighbor emission- 
line free region to account for potential systematic errors brought 
in the continuum subtraction; (3) The equivalent width of broad 
Ha line EW(Ha) > 10A, in order to eliminate potential contami- 
nation from the broad wing of the narrow component. It should be 
mentioned that since we use strict criteria here to select broad line 
AGN, we will correct the incompleteness late through simulations. 

Narrow line objects are classified into SF galaxies, Seyfert 
2 galaxies (Sy2), Low Ionization Nuclear Regions (LINERs) and 
composite galaxies, according to their locations on BPT diagrams. 
We adopt the empirical classification scheme proposed by Kewley 
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EW(Hfi) H4000 

Figure 3. The distribution of H5 absorption line equivalent width (left 
panel) and strength of 4000A break (right panel) for Seyfert 2 galaxies 
(dashed red line), SF galaxies (dashed green line), LINERs (dashed purple 
line), and weak lined objects (solid black line). 



et al. (2006). First, we divide narrow line radio galaxies into star- 
forming galaxies, composite galaxies and 'pure' AGN on [Nn]/Ha 
versus [Om]/H/3 diagram. Second, 'pure' AGN are further sepa- 
rated into 'LINER' or Seyfert according to their locations on the 
[Om]/H/3 versus [Sii]/Ha or [Oi]/Ha diagram (see Figure [2}. For 
137 objects with strong detections in [Om] (with S/N> 10) but 
none in H/3 (S/N< 3), we estimate their H/3 by using the average 
Ha/H/3 of 6.08. Most non-detections are caused by low signal to 
noise of the spectrum around H/3 region or low H/J equivalent width, 
but not due to large extinctions. Their 3 cr upper limits are consis- 
tent with the average value of Ha/H/3 ~ 6.08. 

With the above criteria, we obtain 352 broad lined objects. 
However, among them, 43 objects (40 galaxies, 3 quasars) do not 
have reliable broad H/3 flux (with S/N < 5), 13 ones (11 galaxies, 
2 quasars) have broad line Balmer decrement (Ha/H/3)BL > 10, in 
analogy with type 1.8/1.9 Seyfert galaxies (cf. Dong et al. 2005). 
We reject those objects from further analysis because: (1) the ex- 
tinction correction to the broad or narrow lines are uncertain for 
objects without H/3 flux; (2) the sample becomes very incomplete 
at large Balmer decrements due to their large obscuration. 

Finally, 286 broad lined objects, including 248 from SDSS 
quasar sample and 38 from the main galaxy sample, will be used for 
further analysis. We note that 55 of these Type 1 AGN (46 quasars 
and 9 galaxies) are Narrow-line Seyfert 1 galaxies (NLSl's), ac- 
cording to the criteria of Zhou et al. (2006), with broad component 
of Ha or H/3 > 10 cr confidence level and FWHM(Ha) < 2, 200 km 
s -1 . 35 of them are already included in the NLS1 sample of Zhou 
et al. (2006). 

Among the 3,297 galaxies with detectable emission lines, 
1,771 ones can be un-ambiguously classified, yielding 711 Seyfert 
2 galaxies, 342 LINERS, 216 star-forming galaxies and 502 "com- 
posite" galaxies. The remaining 1,526 weak emission line galaxies 
either do not have sufficient number of measured line ratios to al- 
low a meaningful diagnostic on the BPT diagram (1510 objects), or 
the narrow line decrement Ha/H/3 > 15(16 objects). Most of those 
weak lined objects are likely to be LINERs rather than Seyfert 2 
galaxies because the distributions of their H<5 absorption line equiv- 
alent width and strength of 4000A break are more similar to that 
of LINERs rather than Seyfert galaxies (see Figure[3}. Therefore, 
excluding those weak lined objects will not seriously affect our es- 
timate of type 1 fraction. 

The final sample consists of 71 1 type 2 AGN and 286 type 1 
AGN, selected from original 7810 radio galaxies and quasars with 
z < 0.35 and Pi. 4G hz > 10 23 W Hz" 1 . Their redshift distribution are 
shown in Figure [4] Their locations on the radio luminosity versus 
the narrow Ha luminosity Lu„ are shown on Figure[5] For compar- 
ison, we also add the empirical relation for SF galaxies (Yun et al. 




39 F " i 

20 21 22 23 24 25 26 
log(P..4GHz(W Hz- 1 )) 



Figure 5. Radio luminosity Pi.4ghz versus narrow Ha luminosity L\\ a . The 
points represents for all SDSS selected radio AGN, regardless their radio 
luminosity. The solid line indicates Pi 4GHz _ Lu a relation in SF galaxies. 
Beneath this line, the radio luminosity Pi 4ghz will be barely dominated by 
star formation. The dashed line indicates the 10 23 W Hz~' threshold. 

2001; Kewley et al. 2002). All radio AGN locate on the lower right 
regime, suggesting that there is little radio contamination from the 
star-forming process. This is true, even if the sample is extended to 
a much lower radio power. 

2.4 Black Hole Mass, Nuclear Luminosity and Eddington 
Ratio 

For the Seyfert 2 galaxies, the black hole mass (M.) is estimated 
from the stellar velocity dispersion using the M. - cr, relation ob- 
tained from local host spheroids (Tremaine et al. 2002), includ- 
ing elliptical galaxies and the bulges of disk galaxies. The intrinsic 
scatter of this relation is estimated to be about 0.3 dex. The stellar 
velocity dispersion has been obtained by fitting the SDSS spectrum 
with the IC templates (refer to §2.2). It should be pointed out that 
the stellar light within 3'.' of the SDSS fibre may contain only a 
fraction of galaxy spheroid or have a substantial contribution from 
the galactic disk. In the former case, the correction is usually small 
for most of our galaxies according to the formula of lorgensen et 
al. (1995). In the latter case, the measured cr, may over-estimate 
the stellar velocity dispersion of the bulge component in high in- 
clination disk galaxies, or under-estimate the true value in low in- 
clination disk galaxies. The correction depends on the size of the 
disk and bulge, the bulge disk ratio, the bulge velocity dispersion 
to the circular velocity of the disk, the distance of the galaxy, as 
well as the inclination of the disk. There is no simple formula for 
this correction. Fortunately, most radio selected AGN have large 
black hole masses, thus reside in the galaxies with large bulges. In 
these galaxies, the disk-light contamination within SDSS fibre may 
not be very severe. In addition, the rotational velocity is propor- 
tional to the central stellar velocity dispersion of these large galax- 
ies (Courteau et al. 2008). Therefore, this correction may be less 
severe to our sample. 
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Figure 6. Black hole mass log(M.) versus the absolute optical magnitude 
M r for the radio galaxies. Left panel displays the Seyfert 2 galaxies. The 
orange points represent for the high inclination disk galaxies with axial ra- 
tio b/a < 0.5 and o",/err > 10, the green points for low the inclination disk 
galaxies with b/a > 0.8 and cr./err > 10, and black points for elliptical 
galaxies. The solid line shows the mean M r value in the log(A/.) bins. The 
dotted line show the M r -M. fitting of inactive galaxies suggested by Mclure 
& Dunlop (2002), and the dashed line is referred from Faber- Jackson rela- 
tion of Desroches et al. (2007). The right panel displays the M,--M. relation 
of the absorption line galaxies for comparison. The entire objects in these 
two panels have radio luminosity log(Pi .4GHz/W Hz ) > 23. 



In order to check whether the latter effect is important for 
our sample, we divide the galaxies into the disk and bulge domi- 
nated by using the profile parameters provided by SDSS pipeline. 
We consider it a disk galaxy if the likelihood for the exponential 
profile fit is larger than de Vacouleurs model fit by 0.2 dex, or in- 
verse concentration index petroRSQ / petroR90 > 0.33. Note that 
the concentration index criterion C\ = 0.33 from Shimasaku et 
al. (2001) would induce a 15 ~ 20% contamination from oppo- 
site types of both disk galaxies and elliptical galaxies. But we con- 
sider this contamination acceptable. We separate the disk galaxies 
into high and low inclination groups according to their axial ratios 
(b/a; b/a < 0.5 for high inclination systems, and b/a > 0.8 for 
low inclination systems). Besides, all oblate (b/a < 0.3) objects 
are considered as high inclination disks, regardless of their radial 
profile. Figure [6] shows their distribution on the M r versus M. di- 
agram; where M r is the ^-corrected absolute Petrosian magnitude 
of the host galaxy in the galaxy rest frame by interpolating the five 
SDSS apparent magnitudes with Spline function. For comparison, 
we plot the Seyfert 2 galaxies of our sample in the left panel and 
plot the absorption line dominated radio galaxies in the right panel. 
We find: (1) absorption-line dominated radio galaxies are concen- 
trated on the high luminosity regime, while the spectroscopically 
type 2 AGN spread more to the lower luminosities. (2) M, vs M. is 
tilted for type 2 AGN. It departs significantly from what is defined 
by the absorption line galaxies or the early type galaxies at low M., 
and towards higher luminosity. This can be attributed most likely to 
the disk light contribution. (3) the high inclination galaxies (orange 
points) show systematically lower luminosities than that of low in- 
clination galaxies (green points) at a given M. (i.e. <r,) despite of 
their overlap on the plot, and they locate more closely to that of 
absorption line galaxies. Shao et al. (2007) showed that 'edge-on' 
galaxies are on average 0.8-0.9 fainter in magnitude than that of 
'face-on' galaxies in the r-band due to the dust extinction. This 
seems sufficient to explain the offset between low inclination and 
high inclination disk galaxies observed here without invoking any 
additional effect of a disk component on the a, . 

For type 1 AGN, the virial black hole mass can be estimated 
using the empirical relation between the size of BLR and contin- 
uum luminosity, and the emission line width (e.g., Wandel et al. 
1999; Kaspi et al. 2005; Peterson & Bentz 2006;). We use the broad 
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Figure 7. The [O m] luminosity L[om] versus broad Ha luminosity Z,Ha,Br 
. They are all extinction corrected on the basis of their Balmer decrement 
Ha/H/3. 



Ha luminosity as a surrogate for the continuum luminosity (Wang 
& Zhang 2003), to account for the potential contamination from 
the stellar-light in the weak BLR source and from nonthermal jet 
emission in some cases, using the formula of Greene& Ho (2007): 



M. = (3.0^) x 10 6 



10 42 ergs- 



CD 



iFWHMu 



Mq. 



In section 3.3, this relation will be used also in the simulation of 
selection effect for the potentially overlooked broad emission lines 
in the observed type 2 AGN. The typical uncertainty in this relation 
is estimated to be around 0.5 dex (Vestergaard & Peterson 2006). 

We will use reddening corrected [O m] luminosity as an indi- 
cator for the nuclear power of type 2 AGN (Mulchaey et al. 1994; 
Dietrich et al. 2002; Kauffmann et al. 2003; Haas et al. 2007; Net- 
zer et al. 2006; See Reyes et al. 2008 for an extended discussion) 
for two reasons. First, [O in] emission is less contaminated by star- 
formation process, especially in these massive metal rich galaxies, 
than other lines such as Ha and [On]. Second, Balmer decrement 
can be used as an indicator for the global extinction to the narrow 
emission lines, thus the intrinsic luminosity of [O in] can be recov- 
ered based on the extinction corrected line flux. 

The extinction correction to [Om] luminosity is still a com- 
plex issue. Some authors argued that the some of the Balmer lines 
may be totally blocked according to polarization observation (di 
Serego Alighieri et al. 1997). Because the ionization potential of 
0+ is much higher than that of hydrogen, [O m] emission region 
may be even smaller than that of Balmer lines, and subject more 
severely to the dust extinction. Therefore, extinction derived from 
narrow line Balmer decrements may still underestimate the [O m] 
extinction. However, most of the former arguments are based on 
the directly observed [O in] luminosity rather than on the extinc- 
tion corrected [Om] luminosity. We consider applying correction 
should be better than applying no correction. 

We calibrate the relation between the bolometric luminosity 
and [O in] luminosity with the broad line radio AGN and then ap- 
plied it to the type 2 objects as follows. First, we establish a relation 
between the extinction corrected [Om] luminosity and the broad 
Hey luminosity for the type 1 radio AGN. The extinctions are esti- 
mated based on the Balmer decrements of narrow and broad lines, 
respectively, assuming an intrinsic Ha/H/3 = 3.1. [Om] and broad 
Ho- luminosities are fairly well correlated (see Figure|7](. 

A least-square logarithmic-linear fit yields 

logCLifo.Br) = (1-422 ± 0.069) + (0.984 ± 0.067) x log(L [0 im) (2) 
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Figure 8. The black hole mass M. distribution of type 1 s (estimated from 
reverberation-mapping, blue line) and that of type 2s (estimated from bulge- 
M. correlation, red line). The inset small diagram shows M. distribution 
within restricted L/Z-Edd an d L regime, refer to Figure [121 and §4.2. 
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Figure 9. Eddington ratio log(Z) = log(L/Z.Edd) distribution for type Is 
(blue line) and type 2s (red line). 



, with a scatter in log(LH ff ,Bi) of 0.39. Note that the slope is close 
to one, suggesting a linear relation between the two luminosities, 
which is consistent with Zhang et al. (2008) for radio quiet AGN. 

Second, we estimate the continuum luminosity at 5 100A from 
the broad Ha luminosity using the relation obtained by Greene & 
Ho (2007), for the nuclear light dominated AGN. Finally, bolomet- 
ric luminosity is estimated as L = 9(AL A ) 5lao ^. With the bolometric 
luminosity and the black hole mass, it is straight forward to calcu- 
late the Eddington ratio I = L bo i/L Edd . 

The distributions of the black hole mass for the type 1 and 
type 2 AGN are displayed in Figure[8] On average, the type 1 AGN 
have slightly smaller black hole mass than that of type 2 AGN. 




Figure 10. The extinction corrected [O m] luminosity distribution for type 
Is (blue line) and type 2s (red line). The small diagram shows L[om] distri- 
bution within restricted log(L/Z.Edd) and M. regime, refer to Figure [T2l and 
§4.1. 



The medians are 1.0x10 s Mq and 1.4xl0 8 Mq for the type 1 and 
type 2 AGN, respectively. However, the type 1 AGN are on aver- 
age 1.4 times more luminous than type 2 AGN in the extinction 
corrected [Om] luminosity (Figure IToTl. Jackson & Browne (1990) 
found that radio quasars are a factor 10 more luminous in the [O m] 
line luminosities than that of radio galaxies, and they interpreted 
their results as a larger extinction to narrow line region in radio 
galaxies than that in quasars. The difference found here is much 
smaller because we have corrected the intrinsic reddening using 
Balmer decrement while Jackson & Browne have not. Indeed, the 
average [O m] luminosity of the type 1 AGN is more luminous by a 
factor of 6 than type 2 AGN before the extinction correction. Heck- 
man et al. (1992; see also Meisenheimer 2001) found that the mid- 
to far-IR emission (A, ~ 6 - 50 p.m in AGN rest frame) is 4 times 
stronger on average in quasars than that in NLRG for a 178MHz 
3CR sample (including 42 quasars and 75 NLRG with z > 0.3). 
If the mid-to-far-infrared emissions are isotropic and powered by 
AGN, this implies that some of the NLRG may not be Seyfert 2 
type but be low luminosity LINERs, or the transition of the broad 
line strength occurs in the low Eddington ratio and low luminosity 
region (see §5). 

The Eddington ratio log(£) = log(Lboi/iEdd) distribution was 
displayed in figure [9] We found that the distribution of the type 
1 AGN peaked at log(f) ~ -0.9(€ ~ 0.13) in the range of 
-2 < log(f) < 0, while the distribution of the type 2 objects is 
wider and skewed to a lower log(£) although it is also peaked near 

log(f) 0.9. The more extended tail towards small \og(£) in the 

type 2 AGN may be attributed to a selection effect because the 
broad line components are more difficult to be detected at a lower 
log(£). Note that scatters in the estimation of black hole mass and 
bolometric luminosity will cause a distortion in log(f) distribution. 
Fortunately, the scatter (0.5 dex) of the black hole estimation for 
type 1 using the line width and continuum luminosity is compara- 
ble to the combination of the scatter (0.3 dex) in M. - cr, relation 
and the scatter (0.4 dex) in the bolometric luminosity estimation 
using [O m] luminosity for type 2. Thus, the scatter effect on type 
1 and type 2 distributions is similar. We will discuss this further in 
Section 4. 



3 SELECTION EFFECTS 

Selection effects are introduced during the spectroscopic target se- 
lection and the definition of narrow and broad line AGN sample. In 
this section we will quantify these selection effects as a function of 
nuclear properties for both the type 1 and type 2 AGN. We denote 
the number density of AGN in unit nuclear luminosity (L) and unit 
M. intervals as 0(L, M.). The expected number of AGNs in the nu- 
clear luminosity bin L- L + AL and black hole bin M. - M.+ AM. 
can be written as, 



AN(L,M.) = (f>(L,M.) 



\ p(L s \M.,L)S(M.,L,L s ) (3) 
xdL s V m AL s )AM.AL 



where L s is the luminosity of the band that is used in defining the 
magnitude limit of the sample, e.g., the r-band luminosity of the 
host galaxy dominated targets and i-band luminosity for the nuclear 
dominated targets. p(L"\M., L) is the conditional probability that an 
AGN has U at the given M. and L. Selection function S (M„ L, U) 
is the probability to classify the spectrum as type 1 or type 2; V max 
is the comoving volume corresponding to the maximum redshift 
that an object with a luminosity L s could be detected within the 
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magnitude limit. We have assumed that there is no cosmological 
evolution for both types of AGN within z < 0.35, and we will check 
the assumption in §3. 1 and §4. 1 . Since we are interested in the ratio 
of the two type objects, and the sky coverage for type 1 and type 2 
are the same, we will not consider the sky coverage for the sample 
in the volume calculation. 

Rewriting Eq[3]for the binned data set, we find 



1 h(L, L + AL; M., M. + AM.) 

, :1 V«a,j j£° p(L s \M.,L)S(M.,L,L")dL s 



(4) 



for N AGN in the sample, and h(L, L + AL; M.,M. + AM.) is the 
rectangular function. Due to the additive nature of Eq[4] AGN se- 
lected via exclusive rules can be added up simply. 

With <p(L, M.), we can obtain the density of AGN in the unit 
interval of Eddington ratio i = L\,„\/LEdi and M.\ 

<p(€,M.)= f ^L,M.)d(t-€(L)){de{L)ldL) M .dL (5) 

The density at radio luminosity bin could be wrote as 

HPimbz) = J <P(L,M.)p(P lAGHz \L,M.)dLdM. (6) 

, where p(Piaghz\L, M.) is the conditional probability that an AGN 
has a radio luminosity of Pi.4gh z at a given M. and L. 

3.1 K ra ax 

Different magnitude limits have been used in the selection of dif- 
ferent type spectroscopic targets. The apparent magnitude limits 
for main quasars and FIRST counterparts are ( p8 f < 19.1 on the 
psf magnitude, while for galaxies it is rp ctrosi:m < 17.77 in r band. 
We calculate for each object according to its relevant optical 
magnitude limit and the corresponding optical luminosit)0. Since 
our objects were detected by FIRST, which has a detection limit of 
/i.4GHz > lmjy, a V™* can be estimated from the radio power and 
the radio flux limit for each object. Furthermore, we adopt a cutoff 
in the redshift z < 0.35, so the maximum volume is V0.35. The V max 
is the smallest one among the above three values. 

For objects uniformly distributed in the unverse, the average < 
V/ V max > will be around 0.5. We calculated this value for both type 
1 and type 2 AGN. We find < V/ V max >= 0.55 for the type 1 AGN, 



and < V/V ra 



0.54 for the type 2 AGN. This may be taken as 



an evidence that the number density of AGN increases mildly with 
increase redshift, and we will discuss the evolution effect in §4.1. 



3.2 The Probability Function for Type 2s 

Ignoring the contribution of emission lines to the optical magni- 
tude, the type 2 AGNs are selected according to their host galax- 
ies. In order to quantify the conditional probability p(L'\M.,L) of 
type 2 AGN for a given black hole mass and nuclear luminos- 
ity (or [Om] luminosity), we need to establish a relation among 
the host galaxy magnitude U, the [Om] luminosity Z40111], and 
the black hole mass M.. The strong correlation between the mass 
of the massive black hole and the bulge luminosity of its host 



2 Because we are interested in the ratio of type 1 and type 2 AGNs, rather 
than their comoving density, we do not take the survey area into considera- 
tion. 



galaxy was established for the bulge dominated quiescent galax- 
ies in the local universe (e.g., Hating & Rix 2004) and for the ac- 
tive galaxies and quasars (e.g., Peng et al. 2006). However, there 
is still controversial whether this relation depends on the level of 
nuclear activity (McLure & Dunlop 2002). To check this, we di- 
vide the type 2 AGN into high (log(L[om]/erg s -1 ) > 42) and low 
(log(L[om]/erg s ) < 42) [Om] luminosity groups, and examine 
their distributions on M r - M. diagram. We find that at the same 
M., the [Om] luminous galaxies are brighter than the low Z40111] 
counterparts for only 0.1-0.2 mag on average, which can be ac- 
counted for by their different emission line luminosities. Therefore, 
we will assume that the host's luminosity is independence of the nu- 
clear luminosity at a giving black hole mass. Therefore, we write 
p(L s \M,,L) as p(L s \M.), i.e., independence of the nuclear luminos- 
ity L. 

The galaxies are selected based on their Petrosian magnitudes, 
i.e., U should be the total magnitude of the galaxy. In order to 
quantify p(L'\M.), we need to characterize the relation between the 
black hole mass and the total luminosity of its host galaxy, and the 
scatter of this relation as well. For the known Seyfert 2 galaxies, 
their distributions are already shown in Figure [6] We assume the 
including of the missed Seyfert galaxies does not change this re- 
lation. Because the disk contribution increases as the black hole 
mass decreases, the relation between the galaxy luminosity and 
black hole mass is flatter than the bulge-black hole relation shown 
in Figure(6] For comparison, we also show the M. - M(bulge) rela- 
tion for the inactive galaxies: log(M.) = -0.50(±0.05)M/;(bulge) - 
2.91 (±1.04) (McLure & Dunlop 2002), and the Faber- Jackson re- 
lation of Desroches et al. (2007) for the normal elliptical galaxies. 
Apparently, M r -log(M.) relation of our sample is much flatter than 
the latter relations. To describe quantitatively the distribution of the 
M r as a function of M., we calculate the average value and the 
second momentum of M r over each M. bin, and approximately es- 
timate the p(L s \M.) with a gaussian distribution around their mean 
value. 

We assume the type 2 galaxies are relative easy to identify, 
so S(M.,L,L") = 1 once the object is observed spectroscopically. 
This assumption will fail in two cases: (1) nuclear emission lines 
are too weak, so that some lines used in the spectral classification 
cannot be detected in the spectrum with a typical signal to noise 
ratio of the main galaxy sample; (2) the contamination of the emis- 
sion lines from the HII regions of its host galaxy makes the spectral 
classification as composite or HII type. In the first case, we will 
miss very weak type 2 AGN, i.e., the low Eddington ratio objects. 
But this will not affect our study of the type 1 to type 2 ratios, be- 
cause the detection of a type 1 AGN will require an even higher 
Eddington ratio, and our analysis is limited to the parameter space 
that a substantial fraction of both the type 1 and type 2 can be de- 
tected (see Figure[T2ll. In the second case, we will drop some of the 
type 2 AGN in the late type of galaxies, which tend to have lower 
black hole masses. We will discuss this in §4.2. 



3.3 The Selection Function for Type Is 

The probability of identifying a type 1 object depends on the signal 
to noise ratio of the spectrum and broad line parameters, its pro- 
file and intensity, in a rather complicated manner. To quantify such 
a selection effect, we generate a large number of spectra covering 
the Ha blending and H/3 regimes using Monte-Carlo simulations 
similar to what has been done by Hao et al. (2005) but taking addi- 
tional parameters, like the intrinsic reddening and black hole mass, 
into consideration. These spectra are modeled in exactly the same 
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way as we did for the real data to obtain the emission line param- 
eters, and the type 1 objects are selected with the same criteria as 
described in §2.3. This will give the probability of selecting a type 
1 AGN under different physical parameter regime. 

A simulated emission line spectrum is the sum of three com- 
ponents: the narrow line spectrum (f„a ITOU ,), broad line spectrum 
(f broad) and the noise spectrum (noise). In order to mimic the di- 
verse narrow line spectrum and to avoid complicated noise model, 
we use the observed narrow line spectrum and noise spectrum. In 
order to explore the physical parameters as broad as possible, we 
take the narrow line plus noise spectrum from both the type 1 and 
type 2 AGN. 

The narrow line plus the noise spectrum (/narrow + noise) is 
obtained by subtracting the power-law continuum and Fell mod- 
els for the nuclei dominated objects, or by subtracting the stellar 
continuum model for the host galaxy dominated objects from the 
observed spectrum. It should be noted that this treatment of noise 
spectrum is in-exact for the Seyfert 2 galaxies, because the addition 
of a broad line component and its corresponding nuclear continuum 
would also increase the noise. But its effect is likely to be small be- 
cause a broad line with a height of ten percent of the continuum 
flux would be easily detectable in most spectra. 

For the broad line spectrum (fbmad), we estimated them with 
the parameters of line profile and line flux, which in turn depends 
on the intrinsic broad-line luminosity and the dust extinction to the 
BLR. We denotes the broad line component as, 



fbmad ~ A(X)Kf^ 



(7) 



where A(A) and K are the internal extinction and the scaling factor, 
respectively, and f s bmai is 'standard' broad line spectrum. 

For the type 1 objects, f b s road is the best fitted model for both 
broad H/3 and Ho- after correcting for the internal reddening. For 
the Seyfert 2 objects, we assume that the broad line profile can 
be approximated with a single gaussian, and the line width can be 
estimated from the empirical relation in Eq[TJ with the black hole 
mass estimated from M - cr, relation, and line-flux obtained from 
the observed [O in] luminosity with Eq|2] 

A(A) comes from a set of eight extinction values, which cor- 
responds to a uniformly distributed Hff/H/3 values between 3-10, 
assuming the intrinsic Ho-/H/?=3.1. For a given A(A), a set of ten 
K„ are created randomly following the Ha luminosity distribution 
at the observed [Om] luminosity, i.e., p(Lh„|L[oiii]) (see Figure[7]l. 

To summarize, for each AGN, we build a set of simulated 
spectra that has identical narrow line+noise spectrum as the ob- 
served spectrum but with a series of manually-built broad line com- 
ponents. Assuming the nuclei luminosity scaled with broad line, the 
apparent magnitude and volume limit will change when we vary 
fbmad component of Seyfert galaxies. We retained only those K„ 
and A(A) that makes the optical flux for the simulated spectrum 
within the magnitude limit (r < 17.77 for the host dominated ob- 
jects and i < 19.1 for the nuclear dominated objects). As a result, 
up to 80 spectra are created for each observed AGN. 

We measure the emission line parameters of these simulated 
spectra and determine whether they are broad line AGN as for the 
real spectra. In this way, for each object, we obtain a probability 
as a function of the black hole mass, un-attenuated nuclear lu- 
minosity, the amount of extinction, and apparent magnitude, i.e., 
S(M.,L,Eb-v,L s ). With this selection function, we can estimate 
approximately the distribution of E B _ V from the observed one by 
assuming that the real E B _ V distribution does not depend on the 
black hole mass and optical luminosity. Under this assumption, we 
can correct the observed E B _ V distribution using the above simu- 



■° 4 

I 3 





Ha/Hf? 



Ha/H/5 



Figure 11. Left panel: average correction factor N (Eb-v)/N s ""(Eb-v) 
along Balmer decrement (Ha/HyS)^. Right panel: the observed (dashed 
line) and corrected (solid line) Balmer decrement distribution. 



lated spectrum 

p(E B „ v ) = p obs (E B - V )N'"(E B - V )IN s ""(E B - V ). (8) 

where the p obs (E B - V ) and p(E B - V ) are the observed and real distri- 
bution of E B - V ; N"'(E B - V ) and N"'"(E B - V ) are the number of spec- 
tra identified as Seyfert 1 galaxies and total number of the simulated 
spectra that meet the SDSS targeting criteria. For this correction, 
we have assumed that E B - V distribution is independent of black 
hole mass and nuclear luminosity. The assumption has not been 
fully tested, however, it should not dramatically affect our result as 
far as the correction remains small or the dependence is weak. The 
observed and corrected Ha/H/? is displayed in Figure QT] In com- 
parison with the observed one, which peaks at 3.1, and is identical 
to the value reported by Dong et al. (2008), the corrected Ha/H/J 
distribution shows much more objects at higher extinctions. Inte- 
grating over the distribution of E B _ V , we yield: 



S(M.,L,L S )= S(M.,L,E B - V ,L s )p(E B - V )dE iB - 



(9) 



Because the correction increases fast with the degree of extinction, 
one must keep in mind that the true number may become less reli- 
able at higher extinctions. For this reason, we will limit our analysis 
only to E B _ V < 1.3 (i.e., Ha/H/? < 10), the same criteria we adopted 
in the type 1 AGN selection. 

For quasars, the conditional probability p(L s \M.,L) = S(L S - 
L) and S (M., L, L s ) is independent of L s . For the broad line objects 
selected from the galaxy sample, there is a distribution of L s for a 
given nuclear luminosity and black hole mass due to the significant 
contribution from the host galaxy. Using the relation between the 
host galaxy and the black hole for the type 2 Seyfert galaxies (see 
Figure [6) and the nuclear luminosity from Lh^Bi - L sim relation 
(refer to Greene & Ho 2007), we write 



p(L 5 \M.,L) ■■ 



2ncr 



— f dL"° s 'exp l (LJ 



: '(M.)f 



x exp 



2o-] ) 
(U - L host - Lr(L)) 2 \ 
2cr| ) 



(10) 



where L S, (M.) is the average r-band host luminosity in the 
AGN rest frame for a black hole mass M„ and crj is the scatter 
of this relation; L, (L) is the average r-band nuclear luminosity for 
a given L which is estimated via the broad Ho- luminosity. As we 
already seen, cr 2 is much smaller than cr t , so we can approximate 
the second part of the integrand as a d function. 

We show contours of the average selection function on the 
black hole mass versus nuclear luminosity plane in Figure [L2l The 
selection function strongly depends on the BH mass and nuclear lu- 
minosity [O m] : the incompleteness increases towards larger black 



© 0000 RAS, MNRAS 000, 000-000 



Type 1 Active Galactic Nuclei Fraction 9 




7.0 7.5 8.0 8.5 9.0 9.5 

log(M BH ) 




6.5 7.0 7.5 8.0 8.5 9.0 9.5 

log(M BH ) 



Figure 12. The contour of selection function and the restricted regime of 
log(M./M ) > 6.6,log(L [0 m]/erg s" 1 ) > 40.7 and log(L bo i/L Edd ) > -2.1 
(blue lines). The red circles represent for type 2s, blue circles for type Is, 
purple triangles for LINERs, and green diamonds for composite galaxies. 
The contour levels were set from 0. 1 to 0.6, with interval of 0. 1 . 

hole mass and lower luminosities. At [O m] luminosity below 10 407 
erg s -1 , the selection function is below 0.2, even at black hole mass 
as low as 10 8 Mq. The dearth of type 1 AGN below this limit can 
be well attributed to the selection effects. 



4 THE TYPE 1 FRACTION 

With all the above corrections, we will analyze the fraction of the 
type 1 AGN as a function of the basic parameters, such as black 
hole mass, accretion rate, nuclear luminosity and radio power. As 
we can see in the Figure [T2j due to selection effect, the type 1 
AGN can be detected to a reasonable fraction only in a very lim- 
ited parameter space on either M, versus log(L/L Edd ) plane or 
Aoiii] versus M, plane. In order to avoid potentially large uncer- 
tainties introduced with our correction of selection effects, we will 
limit our analysis to the parameter regimes of log(M./Mo) > 6.6, 
log(L [0 m]/erg s' 1 ) > 40.7 and log(L bo i/L Edd ) > -2.1, that a sub- 
stantial fraction of Type 1 AGN are detected. Using the selection 
function S(M.,L,L S ), V msa and p(L s \M.,L) in the last section, we 
can calculate (f)(M., L) for broad and narrow line AGN. 

4.1 The Dependence of Type 1 Fraction on L {om 

The regime of M. and L[ m] used in this analysis is illustrated in the 
upper panel of Figure[l2]with blue lines. In this regime, log(L[oni]) 
distribution of type Is and type 2s are very similar with means 
log(L [0 M]/erg s" 1 ) of 42.47 and 42.40 for type 1 and type 2 AGN, 
respectively (see inserted small diagram in Figure 10). We inte- 
grate 4>{M., L[oiiij) over either M. in the corresponding range to get 
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Figure 13. The dependence of type 1 fraction on L[oni] for observed den- 
sity (dashed line and diamonds), bias corrected density(solid line and dia- 
monds), and maximum redshift evolution corrected density (green line and 
triangles). The evolution correction is based on extrapolation of quasar lu- 
minosity function of Richards et al. (2005). 

0(Aoni]) f° r both type 1 and type 2 AGN. It is straight forward to 
calculate the type 1 fraction by dividing the bias-corrected type 1 
objects density by the total AGN density in each bin. The results 
are show in Figure [131 

The observed type 1 fraction increases with L[om], which is 
similar to, but with a flatter slope than those given by Simpson 
(2005, hereafter S05) and Hao et al. (2005, hereafter H05) for ra- 
dio quiet AGN. After correcting for selection effects, fx keeps at 
nearly a constant value of 20% over the [Om] luminosity range 
40.7 < log(L[oni]/erg s~') < 43.5. The least x 1 fit for a constant 
fx yields f ~ 20.1% which is acceptable at a probability of 0.92. 
This result is different from S05 and H05, who found that the frac- 
tion of Type 1 AGN increases significantly with the nuclear lu- 
minosity for mainly radio quiet AGN. The difference can not be 
considered (solely) as a difference between radio-loud and radio 
quiet AGN, but be attributed to several different treatments. First, 
S05 and H05 does not corrected [Om] luminosity for the inter- 
nal extinction, while such a correction is included in this work. It 
is known that type 2 Seyfert galaxies show systematically larger 
B aimer decrements of narrow lines than Type 1 Seyfert galaxies 
(Cohen 1983; Gaskell 1984; Rhee & Larkin 2005). This extinc- 
tion would reduce the [O in] luminosities of Seyfert 2 galaxies with 
respect to the Seyfert 1 galaxies systematically. In fact, before ex- 
tinction correction, the type 1 fraction in our sample would also 
increase dramatically from ~ 10% at log(L[ m]/erg s -1 ) ~ 41 to 
~ 90% at log(L [0 iii]/erg s~') ~ 43. Second, Type 2 AGN in the S05 
includes also LINERs and composite type AGN. While it is still 
controversy whether composite type 2 AGN are of similar nature 
as Seyfert galaxies, those spectroscopic LINERs are certainly dif- 
ferent from Seyfert galaxies (Kauffmann et al. 2003; Kewley et al. 
2006; Heckmann et al. 2004). Also we note that most composite 
galaxies are located in relative lower black hole mass and LINERS 
on lower accretion rate regime, while broad lined AGN seldom fall 
in these regions (refer to Figure [T2l. Third, as seen in Figure [12] 
the difference can be attributed at least partly to that we focus our 
analysis to the limited parameter regime. Outside the regime, type 
2 AGN are detected in large amount but almost no type 1 object has 
been detected on the regime with large black hole mass and small 
[O m] luminosity. Finally, S05 did not consider selection effect at all 
while H05 used a different selection function. However, as seen in 
the Figure [12] this changes the [Om] luminosity dependence only 
moderate. This is understandable because we limit our analysis to 
the parameter regime, where the correction is only modest. 

Hitherto we ignored the redshift evolution, but V/V mar test 
does show mild positive evolution with z. We estimate the maximal 
impact of the evolution on our result by using the double power-law 
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Figure 14. The dependence of type 1 fraction on Eddington ratio log(f) = 
\og(L/ Lem) for observed density (dashed line) and bias corrected density 
(solid line). 
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Figure 15. The dependence of type 1 fraction on M. for the bias corrected 
density (/f 01 , solid black line and diamonds), the observed density (f^ bs , 
dashed black line and diamonds), and simulated type 1 density (/j™ 1 , green 
line and triangles, based on the convolved type 2 density). Each bin contains 
the same number of sources. 

form quasar luminosity function 0(L, z) for z < 2.1 derived from 
2dF-SDSS LRG and quasar survey (Richards et al. 2005). Assum- 
ing that type Is and type 2s evolved in the same form, in every {L, z] 
bin, we can apply correction of 0(L, z = 0)/C>(L, z) to both type Is 
and type 2s. The result is shown in figure [T3] (green line and trian- 
gles). The redshift evolution pulls down the type 1 fraction by ~ 6 
percent at high luminosity end, but change little at low luminosity 
end. As expected, this effect does not significantly alter our result. 

4.2 The Dependence of Type 1 Fraction on M. and 
Eddington ratio 

The regimes of M. and t = L bol /L EM used in this analysis are 
shown in the bottom panel of Figure[T2] In this regime, the average 
M. for type Is (~ 10 s 04 Mq) and type 2s (~ 10 s 11 Mq) are indistin- 
guishable (see small diagram in Figure [8j. Within the regime, we 
integrate tp{L, M.) over either L in the corresponding range to get 
4>{M,). Similarly, by integrating of <p(M„ I) (see EqO over M., we 
can obtain ip(C) for both type 1 and type 2 AGN. Note that our anal- 
ysis is restricted to the parameter range illustrated in Figure [T2l blue 
lines), where a substantial fraction of both type 1 and type 2 AGN 
are detected. Thus we can calculate the type 1 fraction straightly by 
dividing the bias-corrected type Is density by the total AGN density 
in each log(^) bin and log(M.) bin. The results are show in Figure 
[JJ]and Figure [731 

From Figure [T4l we find that the type 1 fraction keeps nearly a 
constant value of 25% when Eddington ratio changes from log(^) ^ 
-1.5 to log(^) = 0.0 for the black hole mass in the range of 6.6 < 
log(M./Mo) < 9.5. We also note that the correction to selection 
effects is only modest in whole range of Eddington ratio, and does 
not significantly affect our conclusion as shown in the Figure [T2l 
In fact, before the correction for selection effects, type 1 fraction 
is more consistent with a constant value over -1.5 < log(^) < 0.0. 



At lower Eddington ratios, type 1 fraction becomes very low, and 
the transition occurs at log(^) ~ -2. The transition is likely real, 
rather than caused by under-estimate of the selection function. As 
seen in Figure[l2] most detected objects in this region are massive 
objects, and the value of the selection function is between 0.3 and 
0.5. Lack broad emission line of those objects maybe represent a 
true difference in the BLR properties or over-luminous of [Om] 
in comparison with its true nuclear luminosity in this parameter 
region. 

As mentioned in §2.4, the uncertainty in the black hole es- 
timate will broaden the distribution, thus affect our results about 
the dependence of type 1 fraction on black hole mass. Black hole 
masses for type 2 are estimated using the M. - cr, relation. This re- 
lation has an intrinsic scatter of only 0.3 dex, while the Type 1 BH 
masses are estimated using the Virial estimator, which may affected 
by the inclination of the AGN and accretion rate, and has a typical 
uncertainty of a factor of 0.5 dex. Because this larger scatter, a tail 
to high masses and low mass would be greater for the Type Is than 
the Type 2s when cut-off (such as Z40111]) is applied. When a ratio 
of Type 1 over Type 2 is taken, this would then give an increasing 
Type-1 fraction at high and low black hole masses as observed. 

In order to correct this bias, we convolve the type 2 log(M.) 
distribution with a Gaussian of cr = 0.4, so that type Is and type 
2s have the same scatter. We restrict to objects in jL[om],M.} re- 
gion, as shown in upper panel of Figure Q/2] Based on convolved 
type 2 density p™™ 2 > the simulated type 1 density fraction /j™ 1 = 

Ptypc/(Ptypei + ^type2^ was displayed in Figure[J_5] represented with 
green line and triangles. While the observed type 1 density frac- 
tion /° bs = P° y b p cl /(p° y b p Cl + P°ype2) was represented with dashed line 
and diamonds. Therefore, in every log(M.) bin, we could correct 
the selection effect from different M. measurement uncertainty for 
a factor of / 1 slra // 1 obs . Combining it with selection function and de- 
tection probability, we plot the type 1 fraction /f 01 = (/>i/((/>i + $2) 
dependence on log(M.) in figure [151 with solid line. The least x 1 
fit for a constant f s yields x 2 = 15.54 for 5 degrees of freedom, 
which rules out a constant /i at a confidence level of ~ 99%. From 
the figure, we can see that the type 1 fraction is ~ 30% in higher 
BH mass region (log(M./Mo) > 8), about 10 percent higher than 
that in lower M. region. f\ also rises in the BH mass bin below 
log(M./Mo) < 7.6. This may be due to the loss of type 2 AGN in 
this sample. Low mass black holes are usually hosted in the disk 
galaxies with relative small bugle, such galaxies are more likely to 
have a star forming disk. Due to relative large aperture of SDSS fi- 
bre, SDSS spectrum will encompass also the star forming disk. This 
will shift some type 2 Seyfert galaxies into the composite- type. Be- 
cause we consider only Seyfert 2 type spectra, those type 2 AGN 
will be missed in our Seyfert 2 sample (refer to the "composite 
galaxies" in Figure [T2l represented with green points). In order to 
check the fraction of such objects, we examine the distribution on 
the BPT diagram of type 1 AGN, which are selected based solely on 
the presence of broad lines. We find that only 4.05%(9) type Is with 
log(M./Mo) > 7.6 show composite spectra on BPT diagram, while 
21.3%(16) of broad line AGN with log(M./M ) < 7.6 do. This 
confirms our suspicion that the contamination of star-formation re- 
gion becomes more important at low black hole masses. 



4.3 The Dependence of Type 1 Fraction on Radio Properties 

We plotted the distribution of radio luminosity Piaghz in the left 
panel of FigureQJo] It is clear that more type 2 AGN distribute near 
the lower Piaghz limit 10 23 W Hz~' than type 1 objects, and fewer 
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Figure 16. Left panel: the radio power Pi.4GHz distribution for type Is (blue 
line) and type 2s (red line). The small diagram shows Pi 4ghz distribution 
within restricted log^/L^dd) and log(M.) regime, refer to Figure fL?! Right 
panel: the dependence of type 1 fraction on Pi.4ghz for observed density 
(dashed line) and bias corrected density (solid line). Each bin in right panel 
contains the same number of sources. 



type 2 objects in the high radio luminosity region Pi.4ghz > 10 
W Hz -1 . After limiting to radio AGN in the parameter regime de- 
fined in Figure [T2l the overall distribution is also quite similar. The 
type 1 fraction increases from 15% at log(Pi.4GHz/W Hz" 1 ) = 23 
to nearly 30% at log(P 1 . 4GHz /W Hz" 1 ) = 24 and then flattened to 
log(Pi.4GHz/W Hz" 1 ) = 26 (Figure [16}. Note that the trend is very 
similar before and after correction for selection effects.^ 2 -test rules 
out the possibility of a constant f\ at a confidence level of 99.5%. 
However, when the lowest radio power bin is removed, the trend in 
the rest 5 bins is consistent with a constant f\ at a probability of 
50%. Lawrence et al. (1991) found that type 2 fraction decreases 
with radio power for a lower frequency selected 3CR radio sample, 
which consists of mainly powerful radio sources. Since our sam- 
ple includes only a smaller number of powerful (Pmghz > 10 26 W 
Hz -1 ) radio sources, thus is not adequate to address whether f t rises 
at high radio luminosity. 

Comparing to lower frequency selected samples, our sample is 
more likely affected by radio selection effects. It is general consid- 
ered that radio emission from an AGN consists of two components, 
an isotropic extended lobe component and a beamed jet component 
(e.g., Urry & Padovani 1995). When the AGN is viewed along the 
radio jet, the jet component with a relative flat radio spectrum is 
boosted relativistically, while the lobe emission with a steep spec- 
trum remains the same, i.e., high frequency radio flux is boosted 
along the jet direction more than low frequency radio flux. As such, 
a high frequency survey is more sensitive to the core-dominated or 
compact, flat-spectrum sources than to the lobe-dominated or dif- 
fused, steep-spectrum sources. Because it is general believed that 
the jet emerges along the symmetric axis, the sample will be biased 
to the face-on system, whose optical spectrum is type 1 in unified 
scheme. 

The beaming effect can be checked with radio morphology. A 
radio AGN will appear as a double-lobe source when it is observed 
side-away, and as core-jet or core dominated structure viewed 
along the jet direction due to beaming effect of the relativistic jet. 
Therefore, the core dominance parameter, defined as C = P c /P t 
(P, = P c + Pi), as a rough indicator of the system inclination. If 
radio jet aligns with the symmetric axis of the dusty torus, one 
would expect that type 1 fraction increases with C. However, only 
about 12% sources in our sample have resolved radio structures. 
The number is still too small to allow to reach a firm conclusion on 
this. 



5 SUMMARY AND DISCUSSION 

We have quantified several selection effects on the AGN classi- 
fication and spectral targeting for different combinations of host 
galaxy properties, nuclear properties and dust extinction from spec- 
troscopic samples of radio loud SDSS galaxies and quasars using 
Monte-Carlo simulations. Type 1 AGN have much strong selection 
effect than type 2 AGN. Type 1 radio AGN can be studied only 
in a very limited parameter space on the black hole mass versus 
Eddington ratio or luminosity diagram. After correction for these 
selection effects, we find that in the limited parameter region: (1) 
type 1 fraction is nearly independent of extinction-corrected [O m] 
luminosity and Eddington ratio at log(L[ ni]/erg s~') > 40.7 and 
log(£) > -1.5; at very low Eddington ratio, there are almost no type 

1 AGN, which can not be accounted for solely by selection effects, 

indicating a transition in the broad line strength around log(^) 2 

(see also Hopkins et al. 2009). (2) The type 1 fraction increases with 
the radio power from 10 23 to 10 24 W Hz -1 , then keep flat until 10 26 
W Hz~'. (3) The type 1 fraction is ~ 30% in M. > 10 8 M© region, 
~ 10% higher than that in lower M. region. 

Our result on the luminosity dependence is quite different 
from previous studies for radio loud quasar and galaxies (Lawrence 
1991; Hill et al. 1996; Reyes et al. 2008), or mainly radio quiet 
AGN from SDSS (Simpson 2005, Hao et al. 2005). The difference 
can be mainly attributed to three factors. First, we use an extinc- 
tion corrected [Om] luminosity while previous authors did NOT. 
The average Balmer decrements of narrow components for Seyfert 

2 galaxies (~ 6.08) is significantly higher than for Seyfert 1 galax- 
ies (~ 3.59), which will give arise to a factor of 5 difference in 
[O m] luminosity for two type of Seyfert galaxies. Second, we re- 
stricted our analysis to the parameter regimes, where a substantial 
fraction of objects can be detected according to selection function, 
this reduces the uncertainty caused by the correction of selection 
effects. The final parameter regime used here is similar to the one 
that found by Hopkins et al. (2009) based on a different approach. 
This cuts down the objects with very low Eddington ratios and high 
black hole mass, thus raises the type 1 fraction in the low luminos- 
ity end. Finally, our AGN classification is based on more recent 
work of Kewley et al. (2006), which gives a better separation of 
LINERs and Seyfert galaxies. The new criteria removes quite some 
more LINERs than the previous criteria. Because most of these 
LINERs have lower [Om] luminosity (see Figure [T2l. this raises 
significantly the fraction of type 1 in the low luminosity. 

Although the X-ray observed AGN sample do not suffer from 
above selection effects, however, the fraction of type 1 object based 
on the optical classification into broad and narrow lined objects 
suffers from the aforementioned selection effects. Thus it is proba- 
bly the same reasons caused a luminosity-dependence of obscured 
type-1 (Hasinger et al.2008; Gilli et al. 2007; Fiore et al. 2008). 
Rowan-Robinson et al. (2009) using infrared to X-ray spectral en- 
ergy distribution as a classification for obscured and un-obscured 
AGN, and they reached a similar conclusion as this paper. 

Two very different phonometrical models have been proposed 
for the unification of broad and narrow line radio AGN. The reced- 
ing torus model was first proposed by Lawrence (1991) to explain 
the decrease of the fraction of narrow line objects with the opti- 
cal luminosity in radio loud AGN, and was extended to radio quiet 
objects later (Hill et al. 1996; Simpson 2005; Hao et al. 2005; Sug- 
anuma et al. 2006). The basic idea is that the opening angle of the 
dust torus is larger for more luminous AGN because dust sublima- 
tion radius increases with nuclear luminosity while the height of the 
torus is assumed to be independent of the bolometric luminosity. 
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Thus the broad-line region can be seen over a larger opening angle 
in more luminous objects. On the other hand, Grimes et al. (2004) 
showed that observations are consistent with dual population radio 
sources, 'starved' low luminosity AGN, which do not have a broad 
emission line region, and 'Eddington-tuned' high luminosity AGN, 
without invoking a receding model. In our analysis, we excluded 
these very low Eddington ratio objects 0. Our results show that the 
torus opening angle does not change with the nuclear luminosity, 
thus do not support 'receding torus' model, but are consistent with 
the two population models. 

Very little has been known about the dependence of the type 1 
fraction on black hole mass or accretion rate. Using the ratio of in- 
frared to bolometric luminosity as an indicator of subtending angle 
of torus, Cao (2005) found that the opening angle of torus increases 
with the central black hole for a sample of Palomer-Green quasars, 
but does not correlate with its Eddington ratio. On the other hand, 
Zhou & Wang (2005) argued that the subtending angle of torus de- 
creases with increasing accretion rate based on the equivalent width 
of narrow FeKo' line. It should be noted that the origin of the nar- 
row FeKo' is not clear and the equivalent width is also sensitive 
to the column density of the absorbing material. Our results agree 
with Cao et al. (2005) in the region of log(M./M Q ) > 7.6. 

It is difficult to compare our results with theoretical models of 
dusty torus. Most dynamic-based dusty torus models are not able 
to predict quantitatively how the opening angle of torus changes 
with the black hole mass or Eddington ratio, although some relation 
may be expected. In the dusty cloud model, the physical process 
that maintains a thick torus is still not clear (Krolik & Begelman 
1988; Zier & Biermann 2002). Beckert & Duschl (2004) proposed 
that a geometrically thick torus can be sustained at a high accretion 
rate by the balance of energy dissipation due to cloud-cloud colli- 
sion and heating due to accretion. They found a torus height to ra- 
dius ratio H/R ~ V2GM/( c s M(R)), where M is the accretion rate, 
M(R) the mass within radius R, and c s sound speed. Their model 
has a clear prediction that subtending angle of the torus increases 
with the accretion rate. Our results do not support this. Radiation 
pressure support is suggested by Pier & Krolik (1992), and the in- 
frared radiation and local heating pressure can be very effective in 
supporting a smooth distributed torus. However, an equilibrium so- 
lution can be found only in a relative narrow Eddington ratio range 
(Krolik 2007; Shi & Krolik 2008). It has still to be demonstrated 
that models taken into consideration of more physical processes, 
such as instability, would finally lead to a prediction of a constant 
torus opening angle over a fairly broad Eddington ratios. 

In an alternative model of disk outflow, the entrained dusty 
cold clouds are responsible for the obscuration (e.g., Konigl & 
Kartje 1994). Dopita et al. (1998) suggested that the torus and the 
accretion disk may interacte by accretion-outflow feedback pro- 
cess. Elitzur & Shlosman (2006) argued that a torus is present only 
at nuclear luminosity greater than 10 42 erg s -1 , and its covering 
factor decreases with increase luminosity. However, its argument is 
based on a strong assumption that the size of dusty clouds are scaled 
only with their launching radius. On the other hand, numerical sim- 
ulations showed that luminous and high accretion rate (Lboi/^Edd) 
AGN are likely strongly affected by obscuration in the disk wind 
models (Schurch, Done & Proga 2009). 

It should be noted that there are obscuration sources other than 
dusty torus. The large-scale galactic dust lane and the nuclear star- 



3 Most of low Eddington ratio objects (with L/iEdd < 0.01) are LINERs 
(see Figure [121 . 
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Figure 17. The fraction of high inclination ((b/a) cxp < 0.5) type 2s among 
disk host galaxies. Left panel displayed its dependence on log(M.), and the 
right panel displayed its dependence on log(Pi 4GHz)- Each bin contains the 
same number of sources. 

burst region are certainly responsible for obscuration in some of 
observed type-2 AGN. As we have noticed that a fraction of type-2 
objects in the sample resides in relative edge-on disk galaxies. In 
Fig[T7] we displayed the fraction of high inclination (b/a < 0.5) 
type 2 galaxies among disk host galaxies. The fraction decreases 
quickly as black hole mass and radio luminosity increases. This 
matches almost exactly the variations in f t with M. found in Fig 
[TBI a ~ 10% step function in M. at around 10 8 Mq. Also, the de- 
crease of the edge-on disks fraction in right panel of Fig[l7]may be 
related to the large increase of fi in the lowest radio power bin in 
FigUH These results can be explained by the the presence of a pop- 
ulation of Seyfert galaxies, seen as Type-2 due to galaxy-scale ob- 
scuration. Therefore, the increase of type 1 fraction with black hole 
mass and radio luminosity may be understood as that host galaxies 
moves gradually from disk to early-type galaxies, the extinction by 
the host galaxy decreases. This result offer a very different interpre- 
tation from the rest of the paper discusses. 
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